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The reactions pp/pp — tiH + X represent important channels in the search for the Stan- 
dard Model Higgs boson at the Tevatron and the LHC. In the leading perturbative order 
the cross sections suffer from severe (renormalization and factorization) scale uncertain- 
ties. The next-to-leading order QCD corrections stabilize the cross sections considerably. 
The calculation of these corrections is briefiy reviewed, and a few numerical results are 
discussed. 

1 Introduction 

One of the most pressing open questions in high-energy physics concerns the mecha- 
nism of electroweak symmetry breaking. In the electroweak Standard Model (SM) and 
its extensions, the symmetry breaking is reahzed spontaneously by the Higgs mechanism 
which predicts the existence of Higgs bosons. Thus, the search for Higgs bosons is among 
the first incentives in present high-energy experiments. For the SM Higgs boson, a lower 
experimental mass limit of Mh > 114.4 GeV (95% CL.) has been set by the LEP exper- 
iments [ 0. Moreover, an upper bound of < 196 GeV (95% CL.) [|^ results from a 
fit of the SM parameters in the predictions for electroweak precision data. 

In the near future, the search for Higgs bosons will continue at the pp collider Tevatron 
[ ^ and later at the pp collider LHC [ ^ . Among the various discovery channels for Higgs 
bosons in the intermediate mass range, Higgs-boson radiation off top quarks, pp/pp — > 
ttH -\- X, plays an important role. Although the expected rate is low at the Tevatron, a 
sample of a few but very clean events could be observed for Mh ^ 140 GeV [ |[ . At the 
LHC, ttH production is an important search channel for Mh ^ 125 GeV [^. Moreover, 
analyzing the ttH production rate at the LHC can provide valuable information on the 
top-Higgs Yukawa coupling [ 0| . 

Pure leading-order (LO) predictions for the pp/pp tiH + X cross sections, which 
have been available in the literature [ |^ for a long time, are notoriously imprecise, since 
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Figure 1: A generic set of diagrams (a) for the Born level, (b) one-loop corrections, (c) 
gluon radiation and (d) parton splitting in the subprocesses qq, gg tiH etc. 



they suffer from considerable uncertainties owing to the strong dependence on the renor- 
malization and factorization scales. This uncertainty can only be reduced by including 
higher-order QCD corrections. More recently the complete calculation of the next-to- 
leading order (NLO) QCD corrections to the total cross section at the Tevatron as well 
as the LHC was presented in Ref. [^. These results are in agreement with the parallel 
calculation of Ref. [ |T0| where, however, only the qq annihilation subprocess, which domi- 
nates at the Tevatron, has been taken into account. As expected, in the NLO predictions 
of the cross sections the scale dependence is reduced significantly. 

In this brief article the salient features of the calculation [ ^ ^ of the full NLO 
QCD corrections are summarized. In the next section, the main obstacles in the actual 
calculation and the basic ideas of their solution are explained; a detailed description with 



explicit results is contained in Ref. [|Tl[. Finally, a few examples for phenomenological 
results are discussed, including the scale dependence of the total cross sections at the 
Tevatron and the LHC as well as the transverse-momentum and rapidity distributions of 



the Higgs boson at the Tevatron; more detailed results can be found in Refs. [ P, 11 



2 Features of the NLO calculation 



(i) General remarks 

A survey of typical Feynman diagrams contributing to the processes pp/pp tiH + X 
in LO and NLO is shown in Fig. |1]. In LO the hadron collisions proceed via qq annihilation 
and gg scattering at the parton level, as illustrated in Fig. |l](a). In NLO, virtual correc- 
tions are induced by one- loop diagrams, such as depicted in Fig. |I|(b), and real corrections 
are induced by gluon bremsstrahlung and by parton splittings q ^ qg, q ^ qg in the ini- 
tial state, as shown in Figs. |l|(c) and (d), respectively. Among the one-loop diagrams, the 
pentagons, which have five internal propagators in the loop, are most complicated. The 
idea of the analytical calculation and the strategy for a numerically stable evaluation of 
these pentagons are sketched in the next subsection. The most complicated part in the 
real corrections concerns the extraction of the infrared (IR) singularities and their proper 
cancellation; this is briefly described at the end of this section. 
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Figure 2: Diagrammatical decomposition of the singular part of a pentagon diagram in 
terms of triangle diagrams, where the blobs result from shrinking propagators to a point 
and the fi denote simple kinematical prefactors. 



(ii) Virtual corrections and pentagon diagrams 

The main complications in the calculation of the pentagon diagrams are IR (soft 
and collinear) singularities. The most convenient way to regularize these divergences is 



dimensional regularization. As known for a long time [ |T2[ , in D = A dimensions all 5- 
point functions can be expressed in terms of 4-point functions, simplifying the calculation 
considerably. In order to derive a smiliar reduction of 5-point functions to related 4-point 
functions in D dimensions, first the dimensionally regularized integral E^^'' is translated 
into another regularization scheme that is defined for D = 4. For instance, it is possible 
to endow all massless propagators in the loop with an infinitesimal mass A, defining 
the new integral E^^'^^'^\ which is identical to E^^^ if A = 0. In the next step, a 
related well-defined integral, denoted E^^^^'^\ is determined that possesses the same IR 
singularity structure as E^™^^'^\ This integral is obtained by decomposing the integrand 
of the 5-point function in the collinear /soft limit in terms of 3-point integrands with 
regularization-scheme-independent kinematical prefactors. The explicit construction of 
-^sing^^'^^ is described in Ref. [ |ll]] in detail and illustrated in Fig. ^ for a specific example 
diagrammatically. The difference of the two integrals, E^'^^^^'^'^ — E^^^^'^\ has a uniquely- 
determined integrand and it is regularization-scheme independent, i.e. the limits D ^ 4 
and A — 5> commute in this quantity. In total, we have generated in this way the relation 

E(^) - + 0{D - 4) = ^(--.^=4) _ 4--^'^=^) + 0(A), (2.1) 



where A = on the l.h.s. and D = 4 on the r.h.s. Consequently, solving Eq. ( |2.1| ) for 

E(^\ 

E(^) = E^^^l + [^(-^^^'^=4) - Ei^p''='^^] + . . . , (2.2) 

this integral is expressed (up to irrelevant terms indicated by the ellipsis) in terms of 
four-dimensional integrals and D-dimensional 3-point functions. 

The translation of D-dimensional integrals into four dimensions works not only for 
complete Feynman diagrams, but also for individual scalar and tensor integrals. As shown 
in Ref. [ Q , for D = 4 the scalar 5-point function Eo can be expressed in terms of a linear 
combination of five 4-point functions. Thus, Eq. ( |2.2D expresses Eq^^ in terms of scalar 3- 
and 4-point functions. For the evaluation of 5-point tensor integrals, two entirely different 
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methods have been used. In one calculation the well-known Passarino-Veltman algorithm 



[ is adopted to algebraically reduce the tensor coefficients recursively to tensors of 
lower rank, eventually leading to scalar integrals. This procedure requires to solve a set 
of linear equations for each tensor rank and, in this way, adds a factor of an inverse 
Gram determinant in each step. At the phase-space boundary the Gram determinants 
vanish, since the momenta that span the tensors become linearly dependent. Near the 
phase-space boundary this leads to numerical instabilities that are controlled by careful 
extrapolation out of the safe inner phase-space domains. The second calculation, however, 
avoids the appearance of leading Gram determinants (i.e. Gram deteminants formed by 



four momenta) by using the method of Ref . [ |T4| , which is based on a generalization of 
the strategy [ |12[ for scalar 5-point functions and reduces 5-point tensor coefficients to 
4-point integrals directly. Applying this alternative renders the virtual correction near the 
phase-space boundary much more stable than in the usual Passarino-Veltman approach, 
and the extrapolation from the inner part of phase space turns out to be practically 
unnecessary. The results obtained by the two methods mutually agree with each other. 

(in) Real corrections and IR singularities 

Although the 2 — > 4 scattering matrix element of the real NLO correction o"'"^'^' [see 
Figs. |^(c) and (d)] are quite involved, they still could be calculated by conventional trace 
techniques or by using the Madgraph package [ |15|. The IR (soft and coUinear) singu- 
larities appearing in the phase-space integral of the squared matrix elements have to be 
treated in dimensional regularization and they must be extracted before the numerical 



integration over the four-particle phase space. To this end, a generalization [ |16[ of the 
dipole subtraction formalism [ |l^ to massive quarks has been adopted. In this formalism 
the singularities of the cross section o"''''^' are mapped onto a suitably chosen auxiliary 
cross section a^"^^ so that the difference cr''^^' — o"^"^ can safely be integrated numerically 
in four dimensions. Moreover, a^^^ is still simple enough to allow for an analytical inte- 
gration over the singular regions in phase space. The result of this integration consists 
of LO cross sections dressed by universal functions that contain the singularities. This 
integrated cross section is decomposed into a part af^^ that, defined on configurations 
with LO kinematics, cancels the soft and collinear singularities of the virtual corrections; 
and a second part a^^^ that includes the singularities from initial-state parton splitting, 
which are absorbed in the renormalization of the parton densities. Thus the total NLO 
correction Acr^^° can be written as the sum 



Aa 



NLO 
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_part I ;^-Sub 



(2.3) 



in which each bracket is separately finite. 



3 Some numerical results 



Figure |^ shows the total cross sections for tiH production at the Tevatron and the 
LHC for a fixed Higgs-boson mass of Mh = 120 GeV as function of the renormalization 
and factorization scales, which have been taken at the common value /x. It has been 
checked that no accidental compensations of scale dependences are introduced by this 
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Figure 3: Variation of the LO and NLO cross sections with the renormahzation and 
factorization scales for (a) pp — > ttH + X at Tevatron and (b) pp ttH + X at the LHC 
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Figure 4: (a) Normalized transverse-momentum distribution and (b) rapidity distribution 
of the Higgs boson in LO and NLO for pp ttH + X at the Tevatron (M^ = 120 GeV). 
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identification. In the vicinity of tlie central scale /iq = (2mf + Mj|/)/2, the NLO cross sec- 
tions are remarkably stable with very little variation for fi between ^ fi^/S and ~ 3/io, in 
contrast to the LO approximation for which the production cross section changes by more 
than a factor 2 within the same interval. The relative corrections, which are quantified 
by K factors, K = ct^lo/o'lo (with all quantities calculated consistently in NLO and LO, 
respectively) are nearly constant in the Higgs-boson mass range Mh = 100-250 GeV. 

Although at the Tevatron the cross section is strongly dominated by qq annihilation 
for scales /i ~ /iq, the proper study of the scale dependence requires inclusion of the 
gg, gq, and gq channels. If fj, is chosen too low, large logarithmic corrections spoil the 
convergence of perturbation theory, as indicated by the negative NLO cross section for 
~ /^o/5- Apparently the K factor varies from ~ 0.8 at the central scale /i = /io to 
~ 1.0 at the threshold scale /x = 2/io. As explained in Ref. [ the value of the K factor 
below unity can be understood intuitively in the fragmentation picture proposed in Ref. [ 
although the kinematical limit justifying this approximation is not yet realized at the 



Tevatron. 

At the LHC the major part of the cross section is due to the gg channel, which 
gives rise to increased gluon radiative corrections. For the central scale fiQ the relative 
corrections amount to i^' ~ 1.2, increasing to ~ 1.4 at the threshold value fi = 2/io. 



Again the fragmentation picture [ |T8| is approximately compatible with the complete 
NLO prediction. 

Finally, as illustrating examples for final-state distributions. Fig. ^ shows the trans- 
verse-momentum and rapidity distributions of the Higgs boson for pp tiH + X at the 
Tevatron in LO and NLO approximation. Note that for the distribution in the trans- 
verse Higgs-boson momentum Pt,h, the scale fj, is not taken constant but defined by 
/i^ = Pt// + this choice is more natural for large transverse momenta. The ratio 
of the (normalized) NLO and LO distributions reveals that the mere rescaling of the LO 
distribution by a constant K factor would reproduce the NLO distribution only within 
±10% in the most important ranges. For the Pt,h distribution the increase of the ratio 
d(jNLo/dcLo with increasing pt,h can be traced back to the scale variations of duNLO and 
dcTLo, since fj, rises with increasing pT,i? and daLo decreases with increasing /x faster than 
do"NLO- The respective distributions at the LHC, which are displayed in Ref. [|lll, show 
the same qualitative behaviour, but in the pt,h distribution the ratio daNLO / dcTLo varies 
even more strongly with Pt,h- 



4 Conclusion 



The strong scale dependence of the LO cross sections for pp/pp — > UH+X is drastically 
reduced by the NLO calculation, and the theoretical predictions are stabilized. Thus, the 
NLO cross sections can serve as a solid base for experimental analyses at the Tevatron 
and the LHC. This improvement is also observed for the final-state Higgs transverse- 
momentum and rapidity distributions. 
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